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Abstract 

As shown by Cardy [1], modular invariance of the partition function of a given unitary 
non-singular 2d CFT with left and right central charges cl and cr, implies that the 
density of states in a microcanonical ensemble, at excitations A and A and in the saddle 
point approximation, is po(A, A; cl, cr) = cl exp(27ry / c_lA/6 ) • cr exp(2iry/ crA/6 ). 
In this paper, we extend Cardy's analysis and show that in the saddle point ap- 
proximation and up to contributions which are exponentially suppressed compared 
to the leading Cardy's result, the density of states takes the form p(A, A; cl, cr) = 
/(ciA)/(crA)po(A, A; cl, cr), for a function f(x) which we specify. In particular, we 
show that (i) p(A, A; cl, cr) is the product of contributions of left and right movers 
and hence, to this approximation, the partition function of any modular invariant, non- 
singular unitary 2d CFT is holomorphically factorizable and (ii) p(A, A; cl, cr) / (c^cr) 
is only a function of clA and crA. In addition, treating p(A, A; cl, cr) as the den- 
sity of states of a microcanonical ensemble, we compute the entropy of the system in 
the canonical counterpart and show that the function f(x) is such that the canonical 
entropy, up to exponentially suppressed contributions, is simply given by the Cardy's 
result In p (A, A; c L ,c R ). 
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1 Introduction 

Conformal symmetry as a natural extension of the Poincare symmetry has long been noted 
and extensively studied in the formulation of quantum field theories. In two dimensions, 
the conformal algebra becomes infinite dimensional and admits a central extension c, the 
Virasoro algebra 

[L m , L n ] = (m - n)L m+n + -^n(n 2 - l)S m+n . (1.1) 

As put forward and stressed in the classic work of Belavin-Polyakov-Zamolodchikov [2], 
conformal invariance in two dimensions can be powerful enough to render the conformal 
field theory (CFT) solvable. Being infinite dimensional, unitary irreducible representations 
(irreps) of two-dimensional conformal algebra are also infinite dimensional. Each irrep is 
composed of a highest weight state (corresponding to a primary operator) and its descen- 
dants, obtained by the action of L n , with n < 0, on the highest weight state. A CFT is then 
specified by determining the set of L n 's and its primary fields together with their operator 
product expansions (OPEs). A detailed analysis of 2c? CFTs and references to the original 
papers may be found in [3,4]. 
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In a unitary CFT all the highest weight states and hence all the states of the theory, 
have a non-negative norm. This implies that both the conformal weights A (eigenvalues of 
primary states under L ) of all primary operators and the central charge c are non-negative. 
In this work we will only be concerned with unitary 2d CFTs. 

A classification of 2d CFTs may be arranged based on the value of the central charge c 
and the spectrum of the conformal weights of primary fields. Except for some specific cases, 
such as CFTs with a finite number of primary operators (the minimal models) and WZNW 
models, 2d CFTs are not proven to be solvable [3]. In such cases one would like to know 
if conformal invariance teaches us anything about the spectrum and degeneracy of states of 
the CFT. 

A crude though useful notion is whether the spectrum of the primary operators of the 
CFT is continuous (for singular CFTs) or discrete (for non-singular CFTs) as in the case 
of rational conformal field theories or their subclass of minimal models. String worldsheet 
theory on a D- dimensional flat target space, i.e. a system of D free noncompact bosons, 
is a trivial example of singular CFT. This is due to the fact that as a result of translation 
symmetry of the target space, the states and their conformal weights are also labeled by 
the center of mass momentum which is a continuous parameter. Nonetheless, one may 
compactify the target space on T D and make the spectrum discrete, thus effectively turning 
it into a "non-singular CFT" with a discrete spectrum for primary operators [5]. In the 
pursuit of our investigation, we shall only consider non-singular unitary CFTs or CFTs 
which can be made non-singular by a "regularization" like the example of compactified 
string theory mentioned above. In section 2 we will give a more precise definition of what 
we mean by "non-singular CFT" . 

2d CFTs may be defined on two-dimensional surfaces of various topology, in particular a 
torus T 2 . Although not necessary, in many physically relevant cases one needs to make sure 
that the theory is modular invariant. 1 For the case at hand, where the torus T 2 is specified 
by the complex structure (or modular parameter) r, modular transformations are elements 
of PSL(2, Z) =SL(2, Z)/Z2 which act on r (and simultaneously on f) as 



Modular invariance then demands that Z(t, f) = Z(r', f'), where Z is the partition function 
of the CFT on the torus. In a seminal work [1], Cardy realized and emphasized on the role of 
modular transformations and the restrictions modular invariance imposes on the content of 

1 Modular invariance becomes a necessity if the 2d CFT in question is going to be viewed as worldsheet 
description of a string theory [5]. 
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primary operators of the CFT. 2 In the case of minimal models, modular invariance restricts 
the central charge to be a rational number less than one, explicitly c = 1 — &(j> ~^ , where 
p and q are co-prime integers [3]. Moreover, if we also demand unitarity, the central charge 
is further restricted to c = 1 — ^ q+2 ^ q+3 ^ , where q is a positive integer. The theory of 
free fermions on a T 2 is another example of modular invariant theory, if both periodic and 
anti-periodic boundary conditions are allowed. 

Invariance under modular transformations has been employed in unitary non-singular 
CFTs to specify the density of states at a given conformal weight A. It was shown by 
Cardy [1] that in the saddle point approximation and for A ^> c, the density of states p(A) 
grows exponentially as yfK 




p(A) C ardy ~ exp 2lT\ — . (1.3) 



Cardy's analysis may be repeated for a general 2d CFT with left and right central charges 
cl and cr, at left and right excitations A and A respectively, yielding 



p(A, A) Cardy ~ e 27r V^ . . (1.4) 

One can make two observations from the above Cardy formula: 1) The density of states 
in the saddle point approximation is the product of those in the left and right sectors or, 
equivalently, the partition function is the product of a holomorphic and an anti-holomorphic 
function which may be thought of as partition functions of the left and right sectors, respec- 
tively. In other words, in the Cardy's saddle point approximation, modular invariance implies 
that any unitary non-singular 2d CFT is holomorphically factorizable. 3 2) The density of 
states is only a function of c • A and not of their arbitrary combination. 

In what follows we expand upon Cardy's considerations and by exploiting modular invari- 
ance further we explore the validity of the above two observations beyond Cardy's first-order 
saddle point approximation. In particular, we show that for any modular invariant, unitary 
and non-singular 2d CFT the above two observations remain valid up to exponentially sup- 
pressed contributions, within saddle point approximation. We will discuss the precise meaning 
of these two expressions in sections 2 and 3. We should also mention that the expectation of 



2 PSL(2,Z) can be generated by the action of two independent elements, usually known as S- and T- 
transformations, corresponding to (a,b,c,d) = (0,-1,1,0) and (1,1,0,1), respectively. Our analysis here, 
following that of Cardy's, is based on exploiting S-transformation. In an interesting paper, implications of 
T- transformation for the spectrum of modular invariant, unitary CFTs have also been discussed [6]. 

3 As we discuss in section 5, holomorphic factorizability of a 2c? CFT may also be relevant in the context 
of AdS / CFT and 3d quantum gravity theories [7] . 
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validity of the previous observations is compatible with the results of Carlip [8], where first 
order in perturbation theory beyond the Cardy formula was studied. Our analysis here is 
then an all-orders extension of Carlip's results. 

The rest of this paper is organized as follows. In section 2, we revisit the analysis 
of [8] more closely and extend it to all orders in perturbation theory. In this section, we 
have used the more formal language of path integral. Although a clean method in the 
sense of producing an all-orders result in the saddle point approximation, the path integral 
approach may obscure the physical intuition and some of the assumptions we have made in 
the course of the computation. Therefore, in section 3 we present another method based 
on perturbative, order by order expansion in the saddle point approximation, reproducing 
the results of section 2 in a Taylor series expansion. In section 4, using the expression for 
the microcanonical density of states obtained in sections 2 and 3, we compute the canonical 
partition function of the theory as well as its canonical entropy and compare the latter to the 
microcanonical results. In section 5, besides reviewing our results, we draw our conclusions 
and discuss the implications of our findings if the 2d CFT in consideration is regarded as 
the CFT dual to gravity on an AdS3 background. In two appendices we have collected the 
details of some computations, the results of which have been used in sections 2 and 3. 

2 Extension of Carlip's Formulation 

For a generic (Euclidean) unitary 2d CFT with left and right central charges cl and cr on 
a torus with complex structure r = T\ + 1T2, the partition function is defined as 

Z(r,f) = Tr (e 2lTirLo e" 27rifZ °) = ^ p(A, A)e 2nirA e ~ 2mfK , (2.1) 

A,A=0 

where f = T\ — 2T2, A and A, which we assume to take non-negative values, are respectively 
the spectrum of L and Z , and p(A, A) is the density of states with left and right energies A 
and A. Unitarity of the CFT implies that p is positive definite (negative norm states would 
have contributed to the partition function with negative p). In general, Z(t,t) is neither 
a holomorphic function of r nor holomorphically factorizable, where Z(r, f) factorizes to 
holomorphic and antiholomorphic parts, and the modular invariant partition function Z is 
related to Z as 

Z(t, f) = e-^ iTCL - fCR) Z(r, f ) . (2.2) 

If one dealt with a chiral CFT, or the trace in the partition function were only over 
the chiral sector of the Hilbert space of the theory (i.e. over the Lq = sector), one 
could make use of powerful analytic (more precisely meromorphic) functional properties in 
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combination with modular invariance and express the partition function only in terms of 
its polar part and hence arrive at the Rademacher expansion which completely specifies 
the partition function [9,10]. We note that the Rademacher expansion can also be applied 
to holomorphically factorizable CFTs [11]. However, these techniques are not available in 
general and one may wonder how far one can go relying only on modular invariance, the 
question we explore below. 

In his influential paper [1], Cardy took the first steps into this direction, showing that 
for any matching "high-and-low-temperature" expansion of the partition function one may 
compute the density of states for any given 2d CFT. Cardy formula, among other things, 
implies that around the saddle point the partition function is the product of its holomorphic 
and anti-holomorphic parts. This last statement is equivalent to the fact that the total en- 
tropy of the 2d CFT is the sum of the entropy of left and right sectors. In [8], it was shown 
how one may extend Cardy' s saddle point analysis and compute the leading logarithmic cor- 
rection to the Cardy formula. Remarkably, it was observed how the addition of logarithmic 
corrections to the partition function would not spoil its holomorphic factorizability. In this 
section, we wish to broaden Carlip's method and show that indeed this approximate holo- 
morphic factorizability holds true to any perturbative order around the saddle point in the 
1/r expansion, up to exponentially suppressed corrections. We should stress that in order 
to argue for our case we only resort to (i) modular invariance (ii) unitarity of the theory and 
(iii) that the CFT be non-singular, a term which will be defined momentarily. These are 
indeed very mild assumptions in the sense that they are satisfied by many CFTs. 

We begin our analysis by recalling that the density of states at energies A and A can be 
computed using contour integrals in the two corresponding complex planes of q = e 2mT and 
q = e~ 2nif 

p(A,A) = ^| dqdq^-^Z^q). (2.3) 

Note that the variables q and q are to be treated as independent variables and not as 
complex conjugate of each other. Equation (2.3) is the usual Laplace transform to move from 
canonical to microcanonical ensemble, and thus we regard p(A, A) as the microcanonical 
density of states. Let us also remark that for this Laplace transform to be well-defined, 
or stated otherwise, for (2.3) to follow from (2.1), one should assume that the spectrum 
of the CFT is labeled by a non-negative integer n such that A n — A is a non-negative 
integer, where A is the ground state energy, and similarly for the anti-holomorphic sector. 
This requirement, explicitly having discrete and integer- valued A n — Ao, is what defines a 
non-singular CFT. 4 

4 Note that with this definition any rational CFT [3] is non-singular while the converse is not necessarily 
true. Sec footnote 6 for further comments on this point. 
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Next, we use modular invariance to relate Z(t,t) to Z{— 1/r, — 1/r). Doing so yields 



Z(r,r) = e^^e^^- ^ p(A„,A fl )e - e * . (2.4) 



^n,n=0 



Inserting (2.4) into (2.3) we obtain 

p(A, A) = P( A ^ A «) 7 ( a ' 6 *) 7 ( s > M > ( 2 - 5 ) 



n,n=0 

where 



.ioo(+) /.0(+) j 

/(a A) = - / dre" 2 ™^ 2 ™^ = (-z) / ^ e -W-*r* ? (2 . 6 ) 



o 



with 



= A_ 24' a = A -24' 

-A n + ^, 6 n = -A n + ^. 
24 24 



In our conventions, A n > A m if n > m (and similarly for A's). In particular, the ground 
state energy A is the smallest eigenvalue of L . We stress that having a non-singular CFT 
with a mass gap and discrete spectrum, besides in (2.3), has also been used in arriving at 
(2.5). 

One may perform the integral (2.6) using contour integrals and the result will be inde- 
pendent of the details of the contour by virtue of Cauchy's theorem [9]. It can be shown 
that for a < (a < 0) the integral (2.6) is zero and thus we restrict ourselves to positive a 
(a) only. Depending on the sign of b n , the integral I(a, b n ) is either of the form of a Bessel 
function of the first kind J n (z) or its modified version I n (z), for negative or for positive b n , 
respectively. Thus, using (8.412.2) of [12], we recast (2.6) as 



-2tt 



la b 



s- hiAiry/abn) , b n >0 



I(a,K)={ (2-7) 

^J 1 (47T V ^M), b n <0. 
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Plugging the integrals (2.7) into (2.5), one obtains the recursion formula for p 

I— — A I |£8_A I 
p(A, A) = (2vr) 2 V p(A n , A n ,)^ ^/i(4v™ n ) • I* ^h(4irv n ) 



A A 
^n<- 24 24 



(2tt) 2 V p(A„,A„)^ ^7,(4^)-^ ^(4™ 

' » 9/ ... 11 



L±n> 24 i^n<- 2 4 



24 )^" > 24 



1 Ci 
124 


- A n | 






1 24 


-A n | 






1 24 


- A n | 






1 Cl 
1 24 


- A n | 



- (2vr) 2 V p(A„, A n )^ ^/l(47TO n ) • ^ ^Ji(4tt^) 



+ (27T) 2 V p(A n , A n )^ ^ Jl(47T Mn ) • I* ^Jx^TTVn 

< J It 1) 



1 CR 
1 24 


- A n | 






1 24 


- A„| 






1 Cfl 
1 24 


- A n | 






1 c R 
1 24 


- A„| 



(2.8) 



a a 

24 24 



where 



- A n |(A - g) , ^ = J|| - AnI (A - g) . (2.9) 



Equation (2.8), one of our main results, provides a recursive exact formula for the density 
of states and in this respect it may be viewed as analog of the Rademacher expansion 
for a generic unitary, modular invariant and non-singular 2d CFT which is not necessarily 
holomorphic or holomorphically f actor izable. 

To obtain the expression for the density of states we should solve the above recursive 
equation for p. Recalling the behavior of J\(z) and h{z) for large arguments, 

J^r.J—smfl-z), h(z)~-L=e z , z>l, (2.10) 

V 71 Z \4 / V 27TZ 

one can show that, despite the fact that the last three lines of (2.8) contain J\{z) and are sums 
over infinitely many states with (presumably) exponentially growing weights p(A n , A n ), as 
a result of the oscillatory behavior of J\{z), they are nevertheless exponentially suppressed 
compared to the first line for large u n and v n . The details of this analysis have been gathered 
in Appendix A. Therefore, up to exponentially suppressed contributions, the right-hand side 
of (2.8) is given by its first line, namely 

\£l _ A I \£R _ A I 
P (A, A) = (27T) 2 Yl ^ -U^U n ) ■ ^ ^h^7TV n ) . (2.11) 

i±n< 24 '^ n< 24 

Due to the exponential behavior of I n {z), we can safely deduce that in the saddle point 
approximation A 3> c, A ^> c, only the n, n = terms, corresponding to the maximum 
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value of u n and v n , dominate in the sum (2.11), again up to exponentially suppressed con- 
tributions. 5 Hence in the saddle point approximation we can write the density of states p 
as 



, tt 2 V USD h(S%, 
p(A,A)~(-J Po c L ^-^-c R (2.12) 

/ L Ft 

where po = p(Ao, Ao) is the degeneracy of the ground state which is taken to be equal to 
one, and 

.0 _o„./ c - 



with 6 

~c L = c L - 24A , ~c R = c R - 24A . (2.14) 

We emphasize that (2.12) captures all polynomial corrections to the Cardy formula, 
including its logarithmic corrections. It also establishes what we set to prove in this section, 
namely that any 2d CFT is holomorphically factorizable in the saddle point approximation 
up to exponentially suppressed contributions, i.e. the density of states is the product of the 
density of states of the left and right sectors, respectively. Moreover, as explicitly manifest 
in (2.12), p(A, A) depends on the combination of S Q L and S° R or, rigorously, it is of the form 
p(A, A) = cl ■ cr h(Sl) ■ h(S R ). The prefactor (cl • c R ) may be understood recalling that 
p is the density of states whilst the number of states dj\f = p(A, A)dAdA in the ranges 
(A, A + dA) and (A, A + dA) only depends on combinations of and S R . 

We conclude this section by computing the density of states at energy E. For given states 
labeled by A n and A n (A n , A n > 0), the energy E and angular momentum J are defined as 

A - — — -(E + J) , A- — = -(E-J), (2.15) 
24 2 v ; 24 2 v ; v ; 



5 Recalling (2.10), the ratio of n'th and zeroth terms in the sum (2.11) for large u n or v n is proportional 
to e ~ 47r ("o-Mn) e -47rOo-t;„)^ j n ^ c sac [dle point approximation uo — u n 3> 1 and vo — v n ^> 1. 

6 It is usually assumed that the ground state energy is zero, i.e. Ao = Ao = 0, for which case cl = 
c l, cr — Cr. However, there are interesting and important examples with non-zero Ao. One famous 
example discussed in [8, 13] is the Liouville theory, for which Ao is such that c = 1. Another case with non- 
zero Ao, occurs in the context of AdS3 hairy black holes [14]. The linear dilaton CFT (see chapter 2 of [5]) 
provides another example of such theories. Consider D-dimensional bosonic string on the linear dilaton 
$ = V^X^ background. The central charge of the corresponding worldsheet CFT is c = D + QV^V^ [5]. A 
straightforward and explicit calculation shows that the spectrum of this theory is the same as the one of 
a free bosonic string at = 0, except for the fact that the zero-point energy is shifted by Ao = V^V^/A. 
Nonetheless, as is also stressed in footnote 1 of section 7 of [5], the density of states of this theory is the 
same as the one of a string theory at = 0. This may be understood through our equations noting that in 
this case c = D + 6V% - 24A = D and A - £ = A y=0 + - D+ ^ Vli = A v=0 - 53, which is 
independent. Linear dilaton theory is an example of a non-singular CFT which is not necessarily rational. 
Wc would like to thank the anonymous referee for a comment on this point. 



s 



where 

v 12 J ~ ~ 12 ' 24 ~ 

Next, we note that 

/A x , 1 d 2 Af d 2 Af 

where M is the total number of states. One may then compute dM/dE by integrating 
(2.12) over J. We are only interested in high energy contributions, E ^> ||, ||, for which 
— _E < J < E and hence up to exponentially suppressed contributions we have 

P(E) = ^ = \J E _ E dJ A ) = yV^| 2 I 1 (^yj^sm0)I 1 (2v^^coB0) 



— c to t • 1 ^f Cdrdy ^ + exponentially suppressed contributions 

3 >J Cardy 



(2.17) 



where 



Scardy = ^ \j , Qot = C L + C R . (2.18) 

The details of the computation of the integral (2.17) are given in Appendix B. Thus, it is 
evident from (2.17) that p(E), which is the extension of the standard Cardy formula to all 
orders in perturbation 7 in power oil/E, has the same functional form as the generic density 
of states (2.12). In particular, we note that its dependence on E only appears through the 
c to tE combination. 



3 Cardy Formula to All Orders: Extension of the Saddle Point 
Analysis 

In the previous section, under natural but at the same time strong assumptions of modular 
invariance, unitarity and non-singularity, we computed the density of states for any 2c? CFT 
to all orders in perturbation theory around the saddle point, up to exponentially suppressed 
contributions. In what follows, we employ the saddle point method as used by Carlip [8] 
to go beyond his first-order analysis and reproduce our result of the previous section. In 
subsection 3.1, using the modular S'-transformation f3 — > Air 2 / (3, we first generate Carlip's 
result [8] on the logarithmic correction to the Cardy formula. Then, in subsection 3.2, we 
extend Carlip's analysis to extract the leading-order correction as well as all the subleading 

7 It is an all-orders result in the sense that we are dropping terms of the type exp (— av~E) and cxp(—j3/E), 
with a, (3 > 0. 
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contributions to the Cardy formula in the asymptotic regime where (3E is large. Finally, in 
subsection 3.3, we consider the full partition function, including both energy E and angular 
momentum J. Under the same assumptions as before, we retrieve (2.12). Although the 
calculation is somewhat more involved, we find it helpful for two reasons: It shows more 
explicitly the notion of "up to exponentially suppressed contributions around the saddle 
point" and it provides us with a cross-check of the result derived in section 2. We remark 
that in this section E is the energy of a state of a 2c? CFT on a cylinder whereas A and A 
of the previous section are the conformal weights of states of the CFT on IR 2 and hence they 
differ by st±sa. 

3.1 Logarithmic Corrections 8 

We begin by writing the partition function in the usual fashion as 

Z{fi) = ^exp(-/3£ n ) = I dE p{E) exp(-/3£) , (3.1) 

where the spectral density p(E) is a sum of delta functions with positive integer coeffi- 
cients. In spite of the fact that p(E) is a discrete sum of delta functions, in the asymptotic 
regime where E is very large, the levels are very dense and hence is physically meaningful 
to approximate p(E) with a smooth distribution of eigenvalues such that 

p(E)=exp(a(E)) . (3.2) 

We shall further assume that the error we may incur due to this approximation is signifi- 
cantly negligible compared to the error coming from truncating our expansion in corrections 
to the saddle point approximation at any order of our choice. In other words, we are suggest- 
ing to consider (3.2) as an exact equality for some smooth function cr(E). This last remark 
allows us to rewrite (3.1) as (recall the second equality in (2.1)) 

Z(/3) = J dE exp(a(£) - (3E) , (3.3) 

and consider (3.3) as an equality at any order in the large-/? expansion. 

In order to proceed, we draw on modular invariance to match the low-temperature ex- 
pansion of Z((3) 

Z{/3) = exp(-/3£ ) + «i exp(-£E7i) + 0(exp(-/3E 2 )) , (3.4) 

8 The analysis in this subsection was communicated to the authors by Simeon Hcllcrman through earlier 
email exchanges. 
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as f3 — > oo to its high-temperature expansion 

Z(/3) = exp(-47r 2 £ //3) + fll exp(-4vr 2 E 1 //3) + C(exp(-4vr 2 E 2 //3)) , (3.5) 

as /? — > + . Furthermore, let us make an ansatz for the form of cr(E). In reality, since we 
want to reverse-engineer cr(E) to obtain the correct saddle point expansion, we ought to 
make an ansatz for cr'{E) rather than for u{E) 

a'(E) = E q ~ l ■ (b + hE~ p + b 2 E~ 2p ■■■) , (3.6) 

where q is the leading power of E in a(E) and p is some power greater than zero. Equation 
(3.6) gives either the following high energy expansion for cr(E) 

a(E) = \n(K) + V &m E q ~ mp , (3.7) 

m>0 * r 

if g is not a positive integer multiple of p, or 

<r(£) = \n(K) + b q/p \n{E) + V — — £™ , (3.8) 

m>0 H y 
nv£q/p 

if q is a positive integer multiple of p. In the last two equations, K is meant to be taken as 
a constant independent of E. 

With the stage now set, we are ready to attempt the evaluation of the integral (3.3) in 
the high-temperature expansion (3 — > + in terms of the unknown numbers p, q and b m . The 
saddle point for the integral occurs at the energy E* satisfying 

=/3. (3.9) 

From (3.6) and (3.9) we then obtain the leading-order saddle point equation 

b Q Ef 1 = P, (3.10) 

which yields the following leading-order solution 9 

E* = (/3/b y^. (3.11) 

From (3.11), we learn that q must be less than one since E* must increase as (3 — > + . 
Since at leading order the value of a(E*) — f3E* ought to be equal to ln(Z(/3)), we finally 
have 

1—0 — 1 / V A 4ll 2 En ( 47r 2 (B 1 -E ) \ 

-fco^/T^i + o = -iLpi + o e ^ . (3.12) 



'Since both the right-hand side of (3.10) and E* are positive, so must bo. 
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We draw the reader's attention to the absence in the right-hand side of (3.12) of power- 
law (or logarithmic) contributions in fi. This fact entails that the terms subleading to at 
small (3 vanish to all orders in perturbation theory in fluctuations around the saddle point. 
We shall use this observation to bootstrap the higher-order terms in cr(E). 

As it may readily be seen, the limit /3 — > + in (3.12) predicts that q = | and 



b = 2tt^/^E~ = 2tt y/\E^\, (3.13) 
as Eq = — CL 2^ R i n an Y unitary CFT. Thus, the leading high energy expression for <r(E) is 



which gives the usual Cardy formula for p(E) 



p(E) ^ exp <j 2^/ &±£*!E [ . ( 3 .15) 

In fact, substituting in (3.15) the leading-order saddle point value, we do indeed find 

a(E*)-/3E, = -^!j^, (3.16) 

which shows we have made no mistakes at this order. 

We now move on and compute the first-order correction. We have to be cautious about 
two possible sources of error in (3.12). First, an error might come from including either the 
term -^,E q ~ p = -tz^E^~ p in the expansion for cr(E), if p ^ |, or the term b\ h\{E), if p = |. 
The second possibly relevant source of error is the inclusion of the first correction to the 
saddle point approximation in the expression (3.3) for the partition function Z(j3). Let us 
look at the latter contribution first, which is nothing other than a Gaussian integral over 
fluctuations of E about E*. Defining e = E — E*, and expanding around the saddle point 
e = 0, we find 

iir 2 \E„\ ft 3 , , „ / If- , 



«W*> -PE—L^- + O ^ j , (3.17) 

where 

aieading(^) = 2 b^E 1 ' 2 = A7T^\E \E, (3.18) 

is the leading term in the logarithm of the density of states. Integrating over energies yields 



J dE exp (cri eading (£) - (3E) = exp (^—jp^j * J de ex P {" ig^^ \f + ° (" 



E ( 



o 



exp ('^^1 ) ■ (^ 2 \E \V*/3- 3/2 ) ■ (1 + O(Pfe)). (3.19) 
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Hence, 

Fy/3) -^Cardy -heading— fluctuation -heading— correction— to— a ^(/^/^) ; (3.20) 

where F((3) = lnZ(/3) and 

4ti 2 \E q \ 

-TCardy = ~n ) (3-21) 

heading-fluctuation = In (4tT 3 / 2 | E | 1/2 /3" 3/2 ) , (3.22) 

and -Fieading-corrcction-to-o- is the leading correction to F(j3) due to the inclusion of the b\ term 
in cr(E). In principle, we also ought to have included the contribution to F(f3) coming from 
a shift of E* as a function of /3. However, we shall see that this contribution vanishes at 
leading order. By modular invariance, F((3) must equal hcardy + O (exp(— 4ti 2 (Ei — E ) / (3)) 
as (3 — > + . From this last observation, we learn that the leading term in heading-correction- to-o- 
must precisely cancel -heading-fluctuation in the limit (3 — > + , since the latter goes as In j3, which 
is much larger than the correction of size O (exp(— Ait 2 (Ei — Eq) / (3)) appearing in the right- 
hand side of the equation for modular invariance (3.12). 

Having established that first-order corrections to the saddle point approximation vanish 
and thus do not act as possible sources of error, we turn to computing the leading correction 
to the value of E* when the b\ term is included. The saddle point equation (3.10) is then 
modified to 

iC 1/2 + ^iC 1/2 ~ p = j- , (3-23) 



where bo = 2ffy/\E~o\. In order to find the first-order shift in E*, we separate E* into a 
zeroth-order piece and a correction piece as 

E* = + , (3.24) 

and expand the ^-corrected saddle point equation, treating E* as a small quantity to be 
included only at first order. 10 

Let us rewrite the saddle point equation (3.23) in a more appealing form as 

£*-(l + ^y 2 = S, (3-25) 
more suitable to be expanded as 



bo * J P 2 



in 



This treatment will be justified a posteriori at large (3. 
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With a little algebra on (3.26), we find that the solution for E* is 



E^ = 2^(E^) 1 ~ P + 0((E^) 1 ~ 2P ) . (3.27) 



o 



At this point, we investigate the two possible values p can take which are either p ^ 1/2 
or p = 1/2. There are two corrections to the saddle point value of the exponent a(E*) —(3E*. 
The first correction 

A- (E^y~ P + O ((£fr 2p ) , (3.28) 



2 P 

originates from the inclusion of t^--E* 2 P in the expression for a(E*), whilst the second 

5 — P 

correction 



■ [a' (E^) - (3 + O (a" (E^))} , (3.29) 

is a consequence of the shift in value of E*. However, by virtue of the leading-order saddle 
point equation for Ei°\ the second correction vanishes to the order of interest so that the 
total leading contribution to F(/3) from the b\ term is 

leading — correction — to — a \ V * / \V * / J 

= T -^&i- 2p /^- 1 + O (Z? 4 ^ 1 ) , (3.30) 
2 -p 

for p ~. Thus, we gather that the value of p cannot be smaller than 1/2, otherwise there 
would be a nonvanishing term of order fi 2p ~ 1 in F(/3), which must be absent. Conversely, p 
cannot be greater than 1/2: If that were the case, the ln(/3) term in F(/3) coming from the 
fluctuation integral could not be canceled by the b\ term in cr(E). Hence, it follows that p 
must be exactly equal to 1/2, and that the leading large- E behavior of c(E) must be 

a(E) = Ak^\E \E + 6i ln(£) + ln(if) + 0(E~ 1/2 ) . (3.31) 
The partition function is then 



K f dEexp -P(y/E -bo/0)' 



h 2 



E bl , (3.32) 



where A* is an IR cut-off. It is clear that Z(/3) is independent of A* as we are in the high T 
(low (3) regime: A different choice of A* would only affect the lower part of the spectrum of 
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the theory in consideration. Introducing x = vE — bo//3, (3.32) reduces to 



f°° ( 6n\ 26l + 1 

Z(/3) = K dx2[x + ^\ exp 

26i+l poo 



2Kexp 
2Kexp 



\b 2 ~ 




(bo 


J. 






\b 2 ~ 




(bo 


J. 







26i+l 



0(1) 



(3.33) 



The first term in (3.33), namely expyb^/ (3 J , is just the expression of the entropy as given 
by the Cardy formula. Hence, by matching (3.33) to the /3 — > + limit, we are led to solve 



from which we find that 



bi 
K 



3 

~4' 

2- i / 2 iE r /4 



(3.34) 

(3.35) 
(3.36) 



in agreement with the values found in [8]. 



3.2 Beyond Logarithmic Corrections 

Thus far in this section we have reproduced the logarithmic correction to the Cardy formula. 
However, we wish to go beyond the first-order correction and reproduce the expression (2.17) 
which is valid to all orders in the perturbative expansion in 1/r. We do so by introducing 
yet a new variable y = y/]3x which transforms the partition function (3.32) into 

bj Auv 

Z((3) = 2K^=- [ W dy e~y 2 P{y- (3) , (3.37) 

VOo -'-Auv 

where, recalling (3.8) with p — q — 1/2, 

is a polynomial which can be expanded in ^2 jf y\ uv < as f3 — > + . 
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Because of the Gaussian weight, the energy levels above the cut-off have a vanishing 
contribution to the partition function as Auv — > oo. Thus, we can rewrite (3.37) as 



2K 



fa 

c 3 



Vb 



dy e-y P{y;P) 



J -oo 



2K \lr e ^[ 1+ T, d ^ 



(3.39) 



k=l 



where is a function of bo and b m for m = 2, • • • , k — 1. Modular invariance, as discussed 
in the previous subsection, implies that 

4 = 0, keN. (3.40) 

One can then use (3.40) to determine the coefficients b m for all m, leading to an all-orders 
result. It can be shown that 

bm=^r, (3-41) 

where the constants c m are the coefficients of the asymptotic expansion of the logarithm of 
the Bessel function IAz) 



h{z) 



E 



(-i) fe r(§ + *) 



n^i ^ {2z f fcir (| - k) 

3 3 



cxp 



In (^J2ixz 



21 



\z 16z 2 128z 3 128z 4 



(3.42) 



In order to put our findings in a suggestive form and make contact with (2.17), we rewrite 

(3.43) 



m>2 



as 



ii 



Cardy ) 



Scardy + hi ( — C to t 



7T 



3 ^ 2 m C 

2 111 (^Cardy) ~ ^ ' 



( yCtot 



m=2 

In (S'cardy) + hl( h (Scardy ) ) 



m — 1 V Si 



'Cardy 



m— 1 



or 



?r Ctot 



■/i(S, 



Cardy J j 



(3.44) 



(3.45) 



3 Scardy 

where Scardy and c tot are defined in (2.18). Equation (3.45) is obviously the same as the 
result (2.17) obtained in the previous section . 



L We recall that \E§\ = 



K=h\ ^ and 60 = 2^^^ 
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3.3 Generic Case with E and J 



Ideally, we would like to extend our previous analysis and generalize it to the case in which 
J 7^ and Cl ^ Cr. On general grounds, 



oo OO _ 



a(A,A) = V ^A p ~ mr + V — A g - ns + alnA + /31nA + r/(A,A), (3.46) 
^— ' mr t—' ns 

m=—M n=—N 
m^O n^O 

where 

r/(A,A)= ^ c m , n A™A" s , (3.47) 

with A and A related to E and J as in (2.15). 

According to the lore of modular invariance, we have 12 

Z(t, f) = J dA dA exp (<t(A, A) - r A - f A) = exp j -2tt 2 + j , (3.48) 

where we have assumed A min = E /2 and A min = i?o/2. From our saddle point analysis, we 

are led to set N — M — 1, r = s = |, p = g = 0; also, a 2 ^ = — 27r 2 _E and a 2 x = — 2ir 2 E . 
Equation (3.48) then reads 

Z(r,r) = A dx [°° dx (x + — Y (x + — ) e ("™ 2 -™ 2 ) (3.49) 



r / \ r 



£0 t/ 



r / v r / \ a_i / v a_i 



where Xq = and Xo = — Since we are working with a Gaussian integral, we can 

take both xq, Xq — > — oo. Thus, (3.49) implies that a = f3 = — | since - should cancel on the 
left hand side. 

To go beyond saddle point approximation, we need to show that rj(A, A) = 0. Let us 
consider a term in t] like the following 

c m , n A m / 2 A"/ 2 , m-n^0. (3.50) 

This last term can be written as a function of x and x as 



a-i\ m (_ a_i\™ /a-iY m ( a -i\ n ^ / TX TX I 
c m ,n xH x + — ~ c m , n — P \ >- — > m-n^0 



r / \ r / \ r / \ r / \ a -i a -i 



(3.51) 



2 Note that in this subsection r (r) is equal to —2mr {2ttit) of the previous section. 
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where P is a polynomial and P(0, 0) = 1. After taking the expansion 

exp jc m , n (^-i) (^=i) | ~ 1 + c m>n , (3.52) 
we deduce from (3.48) that 

0, (3.53) 



■m.n 



which implies that the theory is holomorphically factorized around the saddle point. 

One can then readily continue the analysis along the lines worked out in subsection 3.2 
and reproduce (2.12). Since the computation is basically the same as that shown in the 
previous subsection we do not repeat it here. 

4 Canonical vs. Microcanonical Entropy 

In the previous two sections, starting from the canonical partition function we derived the 
expression for the microcanonical density of states through a Laplace transform. Then, 
employing modular invariance of the partition function we fixed the form of the micro- 
canonical density of states, up to exponentially suppressed contributions in the saddle point 
approximation. By means of the usual thermodynamical equations, one can then read the 
microcanonical entropy S m , c , upon taking the logarithm of the density of states. Let us first 
consider the pure imaginary r case, corresponding to the L = L = E sector, for which the 
microcanonical entropy is 

S m . c , = In (p(E)) = In (y Ctot ) + In (^y^) (4.1a) 

= In (^/| Ctot ) + hearty _ \ ln ^Cardy + O ' ( 41b ) 

where Scardy is given in (2.18). We note that (4.1a) gives the microcanonical entropy up to 
exponentially suppressed terms whilst (4.1b) captures only the leading log-correction to the 
Cardy formula. The logarithmic correction was also discussed in [8]. 

Given the density of states p(E) (2.17), one may insert it back into the expression for 
the partition function Z 

Z(f3) = JdE p(E) exp(-(3E) , (4.2) 
and compute the canonical entropy S c 

S c = \nZ-P^. (4.3) 
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The integral (4.2) can be performed using formula (6.620.4) of [12] to obtain 

Z{P) = - ( e ^ T - l) , (4.4) 



2 

where T = is the temperature. One then arrives at 13 

S c = In 2 + y Ctot T + 0(e~^ T ) . (4.5) 

That is, the canonical entropy, in the saddle point approximation and up to exponentially 
suppressed contributions, is completely given by the Cardy formula and in particular there 
are no logarithm or other polynomially suppressed terms, which appeared in the microcanon- 
ical entropy (4.1a). In other words, the functional dependence of the microcanonical density 
of states (2.17), namely ^Ix(S), is such that the Cardy formula is the exact expression for 
the canonical entropy, up to exponentially suppressed contributions. 

The above argument can be readily generalized to the case with non-zero J, for generic A 
and A. As we discussed, the microcanonical density of states (2.12) is the product of density 
of states in the left and right sectors. Inserting (2.12) into (2.1) the canonical partition 
function takes the form 

Z{t,t) = Z l {t)-Z r {t), (4.6) 

where 

Z L {r) = \ (e$'+ _ l) , Z R (f) = \ (e^¥ _ i) . (4.7) 

As the details of the computations exactly parallel those of the previous case, we do not 
repeat them again. The canonical entropy is then given by 

S C = S? + S?, (4.8) 

where 

S c = In Z L - t — , S c =\nZ R - r g _ . (4.9) 

In terms of left and right temperatures Tl = an d T R = the canonical entropy 



(4.8) is obtained to be 
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_2 



S c = 21n2 + y (c L T L + c R T R ) + 0(e-^ T \ e "V^) . (4.10) 



1 One may compute the microcanonical temperature T mc . given the microcanonical entropy: T m \ 

dS„ 



dE 



Using (2.17) we obtain 



7T 



2 



o ^tot-^in.c. 



Scardy-MScardy) 



where Iv{x) is the modified Bcsscl function and Scardy i s given m (2.18). 

14 Recall that for a generic Tl and Tr, the partition function (2.1) may be written as Tr (^e~^ LL ° e~^ R ^ a 
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Again, we see that the Cardy formula gives the exact canonical entropy up to exponentially 
suppressed terms. 

5 Discussion and Outlook 

In this work, we have exploited ^-transformation of the PSL(2, Z) modular group to learn 
more about the density of states of a generic unitary, modular invariant and non-singular 2d 
CFT. Interestingly, we have found that the density of states is the product of those of the left 
and right sectors and, furthermore, it only depends on S — c(A — with c = c — 24A in 
each sector, up to exponentially suppressed contributions in the saddle point approximation. 
Our main result, (2.12), may have diverse interesting physical implications. We discuss some 
of them below. 

1. Noticing that the density of states only depends on So, one can construct a class of 2d 
CFTs with different central charges and different spectra, but with the same So. This 
class of CFTs, for which the densities of states are the same, is relevant for "orbifolded 
CFTs" or CFTs on an R x S' 1 /Zfc orbifold. As an evocative and related argument, 
we observe that any Virasoro algebra at central charge c has infinitely many Virasoro 
subalgebras labeled by integer k at central charge ck. This point has also been noted 
in [16]. To appreciate this last comment, let us consider the algebra (1.1) and the 
subset of its generators L nk 

L n = \L n k, n 0, 

(5.1) 

L o ~ S = I ( L ° ~~ s) ' 

where c = ck. It is then straightforward to see that the set of L n also forms a Virasoro 
algebra at central charge c. The relation between L n and L n could be understood 
through a (non-single valued) conformal map w = z k on the complex plane. This map 
in turn, can be viewed as orbifolding the z-plane by Z^. According to this line of 
reasoning, it can therefore be stated that c(A — ^) = c(A — ^). This is suggesting 
that two conformal field theories with central charge and spectrum (c, A) and (c, A) 

or equivalcntly Tr (e~P H ~ l3fJ - J ^J , where (3 is the inverse of canonical temperature, \i is the chemical potential 
for the angular momentum J, Lo + Lo = 2H, Lo — Lo = 2 J and (3l = 1/Tl, Pr = 1/Tr. Therefore, 

2 _ 1 1 2pL _ 1 1 
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respectively, should have the same density of states up to exponentially suppressed 
contributions. 

In light of the above arguments, it might be possible to relax the requirement of having 
integer-valued spacing of the conformal weights A n (c/. discussions below (2.3)) while 
keeping the discreteness of the spectrum: For theories with (kc,A/k) one might hope 
to find a "dual" CFT with (c, A) (see the example below). In this "dual" picture the 
spacing of the spectrum and the mass gap SM is then 1/k, or in a more suggestive 
form 5M ~ 1/c. 

2. The observation in item 1. may be employed in the black hole microstate identifications 
within string theory. In particular, it may be used to relate black holes within the 
family of (-DO, D6)-branes with different number of D6-branes. A discussion on the 
latter subject may be found in [17]. In this context families of (kc, A/fc) CFTs are 
mapped to each other by T-duality. See also [18] for a related analysis. 

3. Although we mainly focused on the analysis of 2d CFT, one of our main motivations to 
study this problem was the desire to shed new light on the AdSs/CFT2 correspondence 
(see [15] for a review). As first pointed out by Strominger [19], the Bekenstein-Hawking 
entropy of BTZ black holes [20] is correctly reproduced by the Cardy formula (1.4), 
which is nothing but the thermodynamical entropy of a 2d CFT with Brown-Henneaux 
central charge [21] at temperature equal to the Hawking temperature of the BTZ black 
hole. 

As shown in details earlier (see also [8]), Cardy formula should be viewed as the 
first term in the saddle point expansion. The natural question is then: What do the 
corrections to the Cardy formula correspond to in the gravity picture? 

Of course, before employing our results to address questions in AdS 3 gravity, one 
should make sure that the assumptions of modular invariance and non- singularity are 
expected to hold for the CFT proposed to be dual to quantum gravity in AdS3. Natu- 
rally, modular invariance is needed if we want the dual gravity to be compatible with 
toroidal boundary conditions on the Euclidean AdS 3 , see [7, 11]. The non-singularity 
and existence of a mass gap are harder to argue for. For example, for the Dl/Db 
system as discussed in [25], there are some regions in the parameter space where the 
dual theory is expected to have a continuous spectrum above a gap. This brings about 
a pathology in the dual 2d CFT. Nevertheless, in generic points of the parameter space 
of the D1/D5 system the dual CFT is non-singular. Then, it seems plausible to assume 
that the consequences of our findings may be of significance for quantum gravity in 
AdS 3 . This is compatible with the arguments of [7,11]. 
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An idea to tackle the question posed above, as motivated by stringy ct '-corrections 
to supergravity (e.g. see [22] for a review), is that the Bekenstein-Hawking area law 
receives corrections from the higher derivative corrections to the Einstein-Hilbert ac- 
tion [23]. This idea is shown to work for certain four- and five- dimensional black holes 
in higher derivative gravity theories [22]. Accordingly, the corrections to the Cardy 
formula may then play the role of corrections to that action whilst the exponentially 
suppressed contributions may be associated with "stringy or non-perturbative correc- 
tions" which do not admit a semiclassical description. Stated differently, the existence 
of corrections to the Cardy formula implies that the effective semiclassical description 
of a quantum gravity in AdS3 admitting a dual CFT description cannot simply be 
AdS3 Einstein gravity without higher derivative corrections. 15 However, 3d gravity 
has specific features: It was shown by Saida and Soda [26] that any higher derivative 
corrections to Einstein gravity in AdS3, irrespective of the details of such corrections, 
only result in a shift in the AdS3 radius and hence a shift in the Brown- Henneaux cen- 
tral charge. It is obvious that our corrections to the Cardy formula cannot be captured 
by only a shift in the central charge. Therefore, this particular idea does not work and 
our corrections should be understood in a different way. 

In view of the results of [26], the answer to the above question lies within the lines of 
section 4. The Bekenstein-Hawking entropy of a BTZ black hole should be viewed as 
entropy of a system in a canonical ensemble. This is due to the fact that a black hole 
as perceived by an observer at infinity is a system at a given temperature, the Hawking 
temperature. This fact becomes more apparent recalling the Wald entropy formula [23], 
in which entropy is associated with a Noether charge while the temperature is fixed 
to be the Hawking temperature. As we showed in section 4, the Cardy formula (4.5) 
does not receive logarithmic or polynomial corrections in powers of c • T, it is exact up 
to exponentially suppressed contributions. 

4. An interesting and important question put forth in [7], and discussed further in [11], is 
whether there is a well-defined AdS3 pure Einstein quantum gravity, defined by a path 
integral over the AdS3 Einstein-Hilbert action, where the metric is the only dynamical 
field with prescribed boundary conditions [21]. In [11] a careful analysis of the path 
integral for pure AdS3 Einstein gravity is carried out, taking into account contributions 
of all "Brown-Henneaux states" . (Brown-Henneaux states are the boundary excitations 
of AdS3 background which respect the Brown-Henneaux boundary conditions [21]. 
These excitations are localized around the conformal boundary of AdS3 [11].) It was 

15 It has been argued that the log-corrections [8] and in general the whole Rademacher expansion [24], 
may be associated with quantum gravity effects. 
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then argued that this path integral does not have the expected form of the partition 
function of a "physically sensible theory". 

Our findings may have implications on this question/puzzle. If there exists such an 
AdS3 quantum theory (minimal, pure Einstein theory, or otherwise) of gravity and it 
admits a modular invariant, non-singular and unitary dual 2d CFT, then according to 
(2.12) this CFT is holomorphically factorizable and one may relax demanding this last 
feature as a requirement. Note also that (2.12) is obtained from (2.8) by discarding the 
exponentially suppressed contributions in the saddle point approximation and that in 
the derivation of (2.8) no assumption about holomorphic factorizability has been made. 
(Equation (2.8) is an exact result which is true for generic 2d CFTs whose partition 
functions are neither necessarily holomorphic nor holomorphically factorizable.) In 
other words, our equation (2.8) should be compared to equations (5.13) and (5.4) 
of [11], or more precisely to the modular invariant partition function given in equation 
(3.7) therein, which gives the entropy associated with a BTZ black hole computed from 
the partition function of AdS3 gravity (calculated as described above). As another 
outcome of this comparison, we observe that equation (5.13) (or equation (3.7)) turns 
into the partition function of a "physically sensible theory" (cf. discussions of [11]) if 
the coefficients C*a are replaced with p(A n , A n ) of our analysis. The physical meaning 
of this observation and its implications will be analyzed elsewhere. 

5. In this work we mainly focused on the results "up to exponentially suppressed contri- 
butions in the saddle point approximation" . The interpretation of these exponentially 
suppressed terms, especially with regard to the question of AdS3 quantum gravity and 
the (BTZ) black holes entropy for the cases involving holomorphic partition function, 
corresponding to BPS black holes, has been discussed in e.g. [9,27]. Given our analy- 
sis here, it is another interesting question to study these contributions for the case of 
generic non-BPS black holes. 
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A Suppression of the Oscillatory Terms in p(A, A) 



In this Appendix we show that the last three terms in (2.8) are indeed suppressed compared 
to the first term. To this end, we estimate these three terms and compare them against the 
first term. We do so by assuming that our claim holds, that is 



, 2 ~ ~ A(4vtmo) Ii(Attv ) 



p(A, A) ~ (2<kYc l c r ^ ^ • ^ , (A.l) 

u v 



where 



cl(A-|), v = J~c r (A-^), (A.2) 



and then replace the above expression for p into the sums in the last three terms in (2.8). 
Noting the exponential growth of the states by energy, the main contribution to the sums 
comes from the large A states, a region of the spectrum in which the sums could be approx- 
imated by integrals over A. The integrals one needs to compute then take the form 



POO 

l(a,b)= dE hiay/E^ibVE) 
Jo 

dx xIi(ax)Ji(bx) . 



(A.3) 



To compute the above integral, we make use of (8.447.2) and (6.511.1) of [12] 



k=0 

1 

b ~ 



(A.4) 



dx J\ (bx) . 



The integral (A.3) then becomes 

X) 

'fc!(Jfc + l)! V2 



k=0 

where 



^ 1 /ax 2 ^ 1 



/■oo 

c k (b) = dx x 2fc Ji(6x) . (A. 6) 

Jo 



If we have the expression for we can then compute X(a, b). 
From (A. 6), one can then show that 

d 2 d (\ 



m*® = ~ Ck+1 ~ Tb U Cfc(fo) 1 ■ (A - 7) 
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To obtain (A. 7), we have used the following identities for J n {z) 

z—J 1 {z)=zJ {z)-J 1 {z), J x {z) = ~ J (z) . (A.8) 
The solutions to (A. 7) are of the generic form 

c k (b) = ^^d k , (A.9) 

where is a ^-independent parameter to be determined and e is an arbitrary number. 
Inserting (A.9) into (A. 7) we obtain 

4+i = -4(/c + e) (fc + 1 + e)d k , d = 1 . (A. 10) 

We next note that from (A. 6) one can read Ck(b) = b~ 2k ~ 1 ~ 2e Ck{b = 1) if we regulate 
the integral by replacing dx with d 1+2e x. In other words, e should be viewed as a dimen- 
sional regularization parameter which will be taken to zero at the end of the computation. 
Therefore, to leading order in e we have 

d k+1 = (-l) k 2 2k k\{k + l)!di, k > 0, di = -4e, (A.ll) 

and hence 

fc=0 

where we have assumed 6 > a (for a > b the sum is not convergent). Thus, X(a, b) = in the 
e — ?• limit. We would like to comment that from (A. 3) one can show that X(a, b) = ^X(|, 1) 
and, as clearly seen from (A.12), our regularization respects this property. 

We wish to conclude this Appendix by emphasizing that a similar result on the vanishing 
of contributions from the Bessel function Ji(z), was established in Appendix B of [9] using 
a different regularization scheme. 

B Computation of the Bessel Function Integral 

Here, we present a detailed computation of the integral (2.17). Noting that I n (z) = i~ n J n (iz), 
instead of (2.17) one may compute 



J(x,y)= d6J 1 (xsm6)J 1 (ycos9) . (B.l) 
Jo 

Let us sketch out the steps for doing this computation: 
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1. By means of formula (8.535) of [12], we express the J n (z) in (B.l) as 

00 \ k 

J 1 (xsm9) = sm9j2—J k+1 (x) (|) cos 2fe 9 , 

fc=0 

oo , 

Ji (j/ cos 9) = cos E ^ J i+1 (2/) (I) sin 2 ' . 



z=o 



2. Using formula (3.621.5) of [12], we can perform the theta integral 

k\l\ 



d9sm 2l+1 9cos 2k+1 9 



o 2{k + l + l)\ 



to derive 

J^V) = E 2 (fc + Z + l)! J ^ J ^ (I)" (| V 



fc,Z=0 

3. With the help of formula (8.440) of [12] 



■ pir(?-*)^. 



we obtain 



where 



y p,g=l ^ ^ 



P- 1 I' 1 I 1 \fc+Z 



c^ = EE 



( p _ fc -i)!( g -Z-i)!(Ar + / + l)! 

P- 1 I- 1 I i \r+s 

(-i) p+f/ EE 



r!s!(p + q — r — s — 1)\ 

r=0 s=0 \r > 1 l 



4. One can show that 

C„ 



™ (p + g-1)! 
5. Thus, we are able to rewrite (B.6) as 



2 ^ {-lf +q /X\2(p+1) , y y(q+i) 

™J^ (p+l)!(g + l)!(p + g+l)! ^2 J 



2(3 



6. Next, let us consider the expansion of J\(\J x 2 + y 2 ) 

g(,, 2/ ),y?T?j 1 ( v ^T?) = 2^ - tL- _t*-) 

n=0 ' ' 

E(-l)" ^ 1 /X\ 2m /y\2(n-m+l) 

n\ ^ m\(n + l-m)\ \2/ V2/ 

n=0 m=0 v 7 



/ x n 2(p+l) /j,\2(g+l) 

(P + !)!(<? 

p+g>-l 



-2 V 

(p + l)!(g+l)!(p + <z+l)! V27 V27 



7. To compare (B.10) with (B.9), we decompose the above sum into three different regions, 
p, q — 0, ■ • • , 00, p — — 1, q = 0, ■ ■ • ,00 and q = —1, p — 0, • ■ ■ ,00 such that 



Q( X y) = -2T (-) 2(P+1) ( VV{9+1) 

Wx,y) z 2^ (v + ma + 1)Hv + a + 1)l { 2 ) { 



p,q=0 



(p + l)\(q + l)\(p + q + l)\\2J \2 



» (-l)P /X\2( P +1) ~ (-1)9 ^^2( 9 +l) 

+ 2 L ^iw 9 +2 2^ 



(B.ll) 



j p\(p+l)\ \2J ^ V2 

8. Using once more (8.440) of [12], we deduce 

j(x, y ) = _ v^!±Z j^y/gnff) + + 

xj/ y x 

We can now take x, 7/ to be imaginary-valued and finally arrive at 



(B.12) 



I(x,y) = / d01i(xsm 9)I\(y cos 6) = — J(ix, iy) 
Jo 

= ^-±2 h(V^+?) - -h(x) - -h(y) . 
xy y x 



(B.13) 



In our analysis, we are interested in the large x, y limit whereby the last two terms in 
(B.13) are exponentially suppressed compared to the first term and may thus be dropped. 
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